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HEAT TRANSFER AND 'THER MAL STRESSES 
IN SANDWICH PANELS 
By Robert T. Swann 

SUMMARY 


Calculated max imum temperature differences between faces and cal- 
culated ther mal stresses are presented for sandwich panels with a pre- 
scribed linear rate of temperature rise at one face and with the other 
face ins ulat ed - Effects of conduction and radiation are included. Maxi- 
mum temperature differences between top and bottom faces are considerably 
less when both radiation and conduction are considered than when radia- 
tion is neglected. From the calculated data an equation is derived that 
relates the max imum temperature difference to the maximum temperature 
difference when conduction only is considered. An approximate method for 
including effects of radiation in calculations of temperature difference 
is presented. 


INTRODUCTION 


One of the structural configurations being used for supersonic air- 
craft utilizes a sandwich panel for the load-carrying skin. The state 
of stress in such panels can be appreciably altered by thermal stresses 
resulting from aerodynamic heating. In order to calculate thermal 
stresses, the temperature distribution through the panel must be known. 

For solid cores, only heat transfer by conduction need be considered 
and the temperature distribution can be readily calculated by an appro- 
priate Idealization of the sandwich panel into a one -dimensional heat- 
conduction problem. If, however, air spaces are present in the core, 
an appreciable amount of heat can be transferred by radiation, with the 
result that the temperature gradients and the resulting thermal stresses 
are reduced. The heat transferred by radiation can be included In the 
one-dimensional heat-balance equation, but because heat transfer by radia^ 
tion is proportional to the fourth power of the absolute temperature, the 
result is a nonlinear differential equation which is difficult to solve. 

In the present paper equations are derived for calculating the tem- 
perature distribution In any sandwich panel in which conduction is pre- 
dominant. Results of an investigation to determine the effect of 
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including heat transfer by radiation in an analysis of the temperature 
distribution in honeycomb sandwich panels are also presented. The effect 
of radiation is evaluated by comparing the maximum temperature difference 
obtained between the panel faces when radiation is considered with that 
obtained when radiation is neglected, all other factors being unchanged. 
These maximum temperature differences are then empirically related. 

Most of the calculated results were obtained with a Reeves Electronic 
Analog Computer. However, an approximate analytical method for calculating 
temperature when both radiation and conduction must be considered is pre- 
sented in an appendix. Limited experimental evidence of the validity of 
the theory is given. 


SYMBOLS 


A 

AA 

B 

C 


D 

E 

F n,m 

K,m 

G(t) 

g(s) 

h 



area 

solidity of core 
temperature -rise rate of core 
arbitrary constant 
heat capacity of material 
arbitrary constant 
modulus of elasticity 
overall configuration factor 
configuration factor 
arbitrary function of time 
Laplace transform of G(t) 
core height 
diffusivity, k/cp 
thermal conductivity 

Laplace transform • 
inverse Laplace transform 
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P 

Q 

1 

S 

T 

T* 

T 

£ 3 ! 

t 

w 

X 

a 

0 

7 

€ 

0 

A 

t 

P 

G 

T 


perimeter of cells of core 

rate of radiant heat exchange between one point and i 
other points 

rate of heat exchange between two points 
width of square cells in honeycomb panels 
absolute temperature 

absolute temper at ure less initial temperature 
nondimens ional temperature 

difference between the temperatures of the faces of the 
panel 

thickness 

core density 

distance from unheated face of panel 
coefficient of thermal expansion 

roots of the equation 0 tan 0 = ^ ^ 

tp 


constant 

emissivity; strain 
Laplace transform of I 1 
dummy variable of integration 
dummy variable of integration 
density of material 
Stefan-Boltzmann constant; stress 
time 

nondimens ional t ime , 9Kt j h^ 
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Subscripts: 

l, 2,3,4,5 

a,b,c,d,e 

C 

eff 

F 

i 

m, n 
max 
o 


function for configuration factors 

stations as identified in figure 1 

temperatures defined in appendix B 

conduction only- 

effective 

unheated face 

positive integer 

integers between 1 and 5 

maximum 

initial condition 


THEORY 


The problem to he considered is the transfer of heat through a 
sandwich panel, one face of which experiences a prescribed temperature 
history which is a function of time only. For solid-core sandwich panels, 
the only mode of heat transfer which must he considered is conduction. 

If the panel is suitably idealized and suitable boundary conditions are 
selected, an exact solution can be obtained by employing the well-known 
one -dimensional partial differential equation of heat balance which 
governs the transfer of heat by conduction. When air spaces exist in 
the core, significant amounts of heat may be transferred by radiation 
from the hot face to the cooler face and to the core elements. When 
the heat transferred by radiation is incorporated into the differential 
equation for heat balance, the result is a nonlinear partial differential 
equation with variable coefficients, the exact solution of which is dif- 
ficult, if not impossible. The problem can be solved by resorting to an 
approximate analysis which results in a set of nonlinear ordinary dif- 
ferential equations suitable for solution on an analog computer . 

The following assumptions were made in this analysis: 

1. Thermal properties are independent of temperature . 



NACA TO 4511-9 


5 


2. Convective heat transfer is negligible. 

5* The temperature of any plane parallel to the faces is uni form. 

In order to obtain a solution which does not involve unnecessary 
complications , the actual core of the sandwich is idealized as follows. 

It is assumed that the temperatures of the faces of the sandwich are 
respectively uni form; therefore, there will be no transfer of heat by 
conduction in these faces. Because of this, the face material can be 
considered as concentrated heat capacitors located at the ends of the 
core elements. It is also assumed that there is no joint resistance 
between the faces and the core elements, and the thickness of the material 
used to attach the faces to the core is converted on the basis of heat 
capacitance into an equivalent thickness of face material which is then 
included in the concentrated heat capacitors. 

The core is considered as a slab with effective thermal properties 
based on the solidity of the core: 


kef f = k AA 
(cp)eff = cp AA 


K eff = 

cp 


where 

k thermal conductivity of core 

AA solidity of core 

cp product of heat capacity and density 

K diffusivity 


Conduction 

For honeycomb panels having internal surfaces of low emissivity, 
or for panels with solid cores (e.g., foan core or laminated plastic 
core), effects of internal radiation may be negligible. The problem of 
heat transfer in these panels reduces to one -dimensional heat conduction 
in the core, governed by the differential equation 
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1 dT(x,r) = & T(x,t) 


K 3t 




The initial temperature of the sandwich is assumed to be uniform.: 

T(x,0) = T q 

The assumed boundary conditions are 


( 1 ) 


(2a) 


T(h,r) = G(t) 


(2b) 


5t(0,t) _ ^F St(0,t) 
Six KM 5t 


(2c) 


where 

absolute temperature 
initial absolute temperature 

idealized thickness of the unheated face 
time 

distance from the unheated face 
core height 

prescribed function of time 

Equation (2b) follows from the assumption that there is no thermal 
resistance between the core and the face, and therefore the end of the 
core is at the prescribed temperature of the face. 

In order to simplify this investigation, the unheated face was 
assumed to be insulated against heat loss. Equation (2c) is derived in 
appendix A from this assumption, plus the conditions that the face tem- 
perature is uniform and that there is no joint resistance between the 
core and the face. It is also assumed that the face and core are made 
of the same material. In order to extend results based on this equa- 
tion to sandwich panels in which faces and core are made of different 
materials, an effective face thickness should be used. This effective 


T 

To 

tp 

T 

X 

h 

G(t) 
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face thickness is equal to the idealized face thickness multiplied by 
the ratio of cp for the face material to cp for the core material. 

The solution of these equations for the temperature at any point 
is derived in appendix A by means of the Laplace transformation and is 


T(x/r) 


= T„ 


2K y 3 Sin g(l - x) 
h2 sin 23 L 

23 



(5) 


where 3p are the roots of the equation 


3i tan 3i = 


h AA 
tp 


w 


A partial listing of the roots of this equation is given in table I. 
A more extensive table is presented in reference 1. 

For the temperature of the unheated face x = 0, equation ( 3 ) 
reduces to 


I(0)T)=Io + ^g i ^_/ o T a W 

23 


e h d A 


(5) 


If G(t) is a linear function of time, Bt, the following expression 
is obtained by substituting the transform of Bt into the transformed 
equation (A8) before taking the inverse transformation: 


T(0,t) = T 


' 

< T _ h 2 

% J 1 

2h AA \ S2_ 

cos 3 



T K 

w 

h AA 2 

i=L 


) J 



( 6 ) 
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The series in this equation is more rapidly convergent than the one 
that would "be obtained by substituting G(t) = Bt in equation ( 5 ) • 


Conduction and Radiation 

If radiant heat transmission is to be considered, the heat balance 
must include net radiant heat exchange as well as conduction. Net radi- 
ant heat exchange between two isothermal surfaces Aq and A2 (as 
viewed from Aq) can be computed from the expression 

«i ,2 - - Ta 1 *) (T) 


where : 

a Stefan-Boltzmann constant 

e emissivity 

t 

The factor Fq^2 is a function of the geometric relation between sur- 
faces Ai and A 2 . This geometric relation is referred to as the con- 
figuration factor and is defined as the fraction of the total radiant 
flux leaving Aq that is incident on A 2 . (See ref. 2.) For the core, 
the area in the above relation can be expressed as P dx, where P is 
the perimeter of the core at a cross section parallel to the faces; 
therefore equations (l) and ( 7 ) can be combined in the following form to 
include the effects of radiation: 


1 8T(x,t) 
K St 


5^T(x,t) 

dx 2 


creP 


k AA 


t^(x,t ) - r 
^0 


q>q(x,|) tH|,t) d£ - 


q>2(x) T^(h,T } - cpj(x) T^(0 ,t) 


( 8 ) 


The symbols cpq denote configuration factors. The terms in brackets 

represent the net radiant heat exchange at any point; the first term 
gives the heat emitted from the point, the second term is the heat flux 
from all other points of the core to the point In question, and the last 
two terms are the heat flux from the faces to the point. 


The initial temperature (eq. (2a)) and the boundary condition at 
the heated face (eq. (2b)) are unchanged, but the boundary condition at 
the unheated face becomes: 



NACA TN 4349 


9 


,w 


1 ar(0,r) _AA atD(0,r) g£ 
K Sr tp 8 x ktp 


T 4 (0,t) - q^(|) T^U,r) d| - cp^ 0*(h,T) 


( 9 ) 


where cp^ again denotes the configuration factors and the terms in 
brackets represent the net radiant heat flux at the unhe ated face. 

Since no exact methods of solution for equations of this type are 
known, a finite-difference procedure will he used. This is done by- 
taking three stations in the core and one station at each face. (See 
fig. 1.) Now equation ( 8 ) and its boundary conditions can be expanded 
by finite differences to a system of four ordinary differential equa- 
tions. In order to expand the radiation terms of equation ( 8 ) into 
finite-difference form, let F n ^ m be the fraction of radiant flux 
leaving all faces at station n which is incident on all faces at sta- 
tion m, multiplied by the area of all the faces at station n, based on 
unit area of the face plates. Expansion of the radiation terms in equa- 
tions ( 8 ) and ( 9 ) into finite differences then gives for each element n 
a term 


Q n = -ere 

m=l 




(10) 


Therefore, in terras of finite differences, equation ( 8 ) and its boundary- 
conditions become 


T 1 


= T 0 



dT P 

— - = 2 T l - 3T2 + Tj 
dr 


ieh 


AA 






dT 


— ± = To - 2Ti + Tk 

dr 2 5 4 3h AA m: 


L_ *5,*(vs k 



(Ha) 

(lib) 

(He) 
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3T 4 + 2T«j - 


geh 
3 k AA 



m=l 



(lid) 


3 T 5 

dr 




h AA 
3tF 





(lie) 


where 


rr _ 9Kr 


T = 


and the initial conditions are 


t = 0 T n = T 0 

Equations of this type are readily solved on an analog computer equipped 
to generate fourth powers of the dependent variable. 

The solution of equations (ll) can be made considerably more gen- 
eral by using a nondimens ional temperature parameter T n which is 
defined by the following relation: 


T n = 7 


fak AA 

V rch , 



( 12 ) 


where 7 is a constant which may be selected so that the problem falls 
in the desired computer range of temperature. 

Substituting this expression for T n into equations (ll) and 


dividing by 



gives 


Ti = T 0 + 



7 3 kAA 
V creh , 


( 13 a) 
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dlo - — — 

= 2D X - 3^2 + T3 



= T > 2 - 2T 5 + Ti,. 


<33?k _ _ _ 

— - = T 5 - 3T4 + 2T 5 
<1t 


dT 5 

d.T 


h AA 
3tp 




- r 3 F 2 , m (i2 4 - <C) 

m=l ' 

- r 3 Si ^ ) 

m=l 

- r 3 *4,m( 5 4 4 - 

m=l 

- r 3 ^ „,fe> - tM 


with the initial condition 


( Tn )t=o 



(15b) 


(13c) 


(13d) 


(I3e) 


An approximate procedure for calculating the temperature of the unheated 
face is presented in appendix B. 


Thermal Stresses 

Temperature differences "between faces of sandwich panels cause 
unequal expansion of the faces and result in thermal stresses. An indi- 
cation of the magnitude of these stresses can readily be obtained if the 
following assumptions are made: 

(1) Material properties are independent of temperature. 

(2) The core is rigid in shear. 

(3) The faces remain plane and all edges are free to expand. 

(4) The core takes no load. 
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The stress- strain relations are 


°i 

€ - gt + an 


C^TTI 

= — + aHV 

% 


where 

e strain 

cf-^ stress in heated face 

crp stress in unh eated face 

E modulus of elasticity 

a coefficient of thermal expansion 

Tp temperature of the unheated ■ f ace, Ttj or T(0 ,t) 

For equilibrium, 


t]_0-j_ t tpOjp — ^ 

Therefore, the stresses in the faces are 


a n = - 


Ea AT 


°F 


t x 
1 + — 
tjp 


Ea AT 


1 + 


tp 

*1 


where : 


( 14 ) 


( 15 ) 


(16) 


AT = T]_ - Tp 
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RESUITS and discussion 


Solutions were obtained for all co mb inations of the following 
parameters : 


h, in 

AA 

tp, in 

B, °R/sec 

cp, Btu/(ft5)(°R) . . . . 

ip op 

k°'(Btu)(ft)/(hr)(°R)(ft2) 


. . 0.2, 0.3, 0.4, 0.5 
0.02, 0 . 025 , 0 . 03 , 0.04 
0.01, 0.02, 0 . 03 , 0.05 

.... 5, 10, 20, 50 

50 

500 

12.5 


In all cases the face thicknesses listed here were modified in the cal- 
culations to include 0.002 inch of face material as the thermal equiv- 
alent of the braze or other material used to fasten the core to the face. 
For cases in which radiation is considered, the core is assumed to con- 
sist of l/4-inch- square cells with eml.ssivity equal to 0.8. Configura- 
tion factors used in the calculations are given in table II. 


Typical results are shown in figure 2. The curve identified as Tp 

is the prescribed temperature of the heated face, which in this case 
experiences a temperature rise B of 10° R per second. The curve 
marked Tp q is the temperature response of the unheated face calculated 

from equation ( 5 ), which applies when radiation is neglected. A com- 
parison of this curve with results calculated numerically and by analog 
computation from equations ( 13 ) without the radiation terms revealed 
differences too small to be shown on this plot. The curve identified 
as Tjp is the temperature of the unheated face calculated from 
equations ( 13 ) • 


As would be expected, the temperature difference between the heated 
face and the unheated face is more when radiation is neglected. The tem- 
perature difference is essentially constant at its maximum value after 
150 seconds when radiation is neglected. (From eq. (6) it is obvious 
that the maximum occurs as time approaches infinity. ) 


The maximum temperature difference with conduction only can be cal- 
culated by subtracting equation (6) from equation (2b). As time becomes 
large, this difference becomes 



(IT) 
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The maximum temperature difference with radiation and conduction was cal- 
culated on an analog computer from equation (ll). When these maximum 
temperature differences were plotted with (ATmax)c as abscissa and 

A T-mfly as ordinate for each case listed at the beginning of this section, 
all the points lay within the envelope shown in figure 3 • 

Since all the points lie in a narrow band, regardless of the values 
•of the parameters, it appears that the effect of radiation depends chiefly 
on the magnitude of the temperature difference which would be obtained 
by considering conduction only. 

The dashed line in figure 3 represents the data in the interval 

120 < (^Tinax)c < 3,000 

with a maximum error of about 10 percent. 

An approximate equation for this line is 



If equation ( 17 ) is substituted for (ATmax)c in equation (l8), the 

following expression is obtained for the maximum temperature difference 
when both radiation and conduction are considered: 



Dimensionless Analysis 

Calculated maximum temperature differences are presented in dimen- 
sionless form in figure 4. The max imum dimensionless temperature dif- 
ference is plotted as a function of the ratio h AA/tp for several 
values of dimensionless temperature-rise rate.' Each plot is for constant 
initial dimensionless temperature, since 



Also, each plot is for a fixed value of the ratio h/S, where S is 
the width of the square cell. 
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On the logarithmic scale used, the curve of maximum dimensionless 
temperature difference varies linearly with the ratio h AA/tp, and the 

slopes generally decrease as the temperature-rise rate increases. How- 
ever, it was not possible to derive a mathematical relation between the 
slopes of the t emp erature-difference curves and the temperature-rise 
rates which would lead to a simpler empirical approximation to the maxi- 
mum temperature difference than that given by equation (l9). Figure 4 
can be used to determine maximum temperature differences with greater 
accuracy than is possible from equation ( 19 )* 


Correlation of Theory and Experiment 

The theory from which the results given in the present paper were 
obtained assumes that no heat is lost from the unheated face of the 
sandwich. There is a scarcity of data on sandwich panels tested under 
conditions which match the theoretical assumptions. However, in a test 
made in the Langley Structures research laboratory a light-weight steel 
honeycomb sandwich panel experienced a temperature-rise rate of 
19 . 2 ° R per second at the heated face and experienced a radiant heat 
loss at the unheated face of the sandwich. Heat was transferred from 
the unheated face of the sandwich by radiation across an air gap to a 
steel backing plate. The sandwich panel properties were tp = 0.005 inch, 

w = 7.2 lb/ft^, and h = 0.3 inch. Test results are shown by the data 
points in figure 5* The prescribed temperature of the heated face is 
shown by the curve identified as Tj_. In order to compute the tempera- 
ture histories of the sandwich faces, it was necessary to modify the 
analysis to account for the radiant heat loss at the unheated face. It 
was therefore assumed that the unheated face of the sandwich panel radi- 
ated to a plate which remained at room temperature. The emissivities 
of the backing plate and the unheated face of the sandwich were assumed 
to be 0.6 and the internal emissivity of the sandwich was assumed to be 
0.8. The calculated temperatures of the unheated face of the sandwich 
are shown by the solid line identified as Ty and are in excellent 
agreement with the test results for the first 50 seconds. After that 
time, the temperature of the backing plate used in the test began to 
rise and thus invalidated the theoretical assumption of constant tem- 
perature and permitted no further basis for comparison. In order to 
show the effect of radiation, a curve of calculated temperature is 
included for the case in which radiation is neglected. 


Thermal Stresses 

In order to give more meaning to the magnitudes of the temperature 
differences between the faces of a sandwich panel, a simplified elementary 
stress analysis was made. If it is assumed that the panel remains flat. 
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that the core takes no load, and that there is no external restraint to 
expansion, the thermal stresses for equal face thicknesses are, from 
equations (l 6), 


Ea A3? 


with the heated face in compression and the unheated face in tension. 

Of course, the analysis overestimates the thermal stresses which actually 
exist; nevertheless, it provides a convenient means of evaluating the 
effect of internal radiation. For steel the product Ea is about 500 
and thermal stresses are 



a = ±150 AT 


( 22 ) 


Therefore, a plot with maximum temperature difference as the ordinate can 
be readily converted into a plot with maximum thermal stress as the ordi- 
nate. If the temperature differences in figure 5 are converted into 
thermal stresses, the lowest stress obtained is 15,000 psi. This result 
indicates that large thermal stresses are to be expected when sandwich 
panels are heated. It appears from equation (19) that high core density, 
thin unheated faces, and thin panels will be necessary in order to prevent 
these stresses from becoming prohibitive. 


However, if the ratio of face thicknesses tq/tp is increased, 
compressive stresses diminish and tensile stresses increase. Compres- 
sive stresses can be reduced by any desired amount without increasing 
tensile stresses by more than a factor of 2. This effect will now be 
considered for a specific example. Take a steel honeycomb panel with 
h = 0.5 inch, B = 20° F/sec, and w = 20 lb/ft5. Assume that the panel 
is to be used in an application requiring a total face thickness of 
0.1 inch; that is, tq + tp = 0.1. The temperature information neces- 
sary for evaluation of the thermal stresses in this example can be cal- 
culated from equation ( 19 ) - It should be noted that, since the temper- 
ature of the heated face is prescribed, its thickness does not affect 
the temperature difference. 3?he face thicknesses given here are the 
actual face thicknesses and do not include a modification to account 
for braze material. The stresses are given by the equations 


0 l,max = “5.j000 AT^x tp 


3 F,max 



* 


= 5,000' ATjaaxCO.l 


( 25 ) 
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Thermal stresses calculated from these equations are plotted against 
unheated-face thickness in figure 6 . If the two faces have the same 
thickness, the thermal stress in each face is 75,000 psi. Decreasing 
the unheated-face thickness from 0.05 inch to 6.636 inch and increasing 
the heated-face thickness to 0.064 inch reduces the stresses in the 
heated face from 75,000 psi to 46,000 psi. Stresses in the unheated 
face increase from 75,000 to 80,000 psi. Further reduction of the 
unheated-face thickness will reduce thermal stresses in both faces. 


CONCLUDING REMARKS 


Calculations of temperature distributions in sandwich panels, even 
when internal radiation is present, are not difficult if analog computing 
equipment is available. If internal radiation is present, the accuracy 
of an analysis of maximum temperature difference between the top and 
bottom faces of the panel, ' based on conduction only, decreases as the 
temperature difference increases. 

If the temperature-rise rate of the heated face of the sandwich is 
linear, a simple empirical relation exists between the maximum tempera- 
ture difference between the top and bottom faces of the panel and the 
maximum temperature difference calculated when only conduction is con- 
sidered. Approximate temperature differences can be calculated by con- 
sidering the separate effects of radiation and conduction. 

Even with relatively moderate heating, large temperature gradients 
exist in sandwich panels and cause large thermal stresses. For a given 
total face thickness, however, the magnitude of thermal stresses can be 
materially reduced by decreasing the unheated-face thickness while 
increasing the heated-face thickness by an equal amount. 


Langley Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Langley Field, Va., June 9, 1958* 



18 


NACA IN k3k9 


APPENDIX A 

CONDUCTION OF HEAT IN SANDWICH PANEIS 


The differential equation for the linear conduction of heat is 


1 8t(x,t 
K dr 


1 _ 8 2 T(x t t ) 
Sx 2 


(Al) 


The initial condition is T(x,0) = T 0 . The temperature of the heated 
face, x = h, Is prescribed as: 


T(h,r) = Ct(t) 


(A2) 


The rate at which heat leaves the core and enters the unheated face Is 


q - k AA 3F19 j-T1 (A3) 

CK 

Inasmuch as It is assumed that the external surface of the unheated face 
is insulated, the temperature-rise rate of this face is proportional to 
the rate at which heat leaves the core: 


,q = cptp 

Equating these expressions for q 


ar(o,T) 

8t 

gives 


8T(0,r) = t F 8t( O.t) 
Sx K AA 3t 


(a4) 


(A5) 


These equations are readily solved by means of the Laplace transforma- 
tion. Let 


T'(x,t) = T(x,t) - T 0 
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The transformed equations are 


d2e ( x ^ - 1 0(x,s) = 0 (A6a) 

dx2 K 


8(h,s) = g(s) 


(A6b) 


59 C ° - ; -gl = fin 0(0, s) (A6c) 

dx KAA 

where 

9(x,s) = Ljr'(x,T)J 
g(s) = L {g( t ) J 

The solution of equation (A6a) is 


0(x,s) 



(AT) 


where C apr i D are constants to be determined from, equations (A6b) 
and (A6c) . When C and D are evaluated, the transformed equation is 



An inverse transform of the equation can be expressed immediately as 
a convolution of the arbitrary heating function with the other terms 
of the equation: 
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T'(x,r) = 



(x < h) (A9) 


where L"-'- indicates the inverse transform. The expression within the 
brackets is single valued with simple poles, and by a formal application 
of the inversion integral it is found to be 


_Kf3ifr 

2K y~ P sin p(l - x) c h2 
h 2 i=I sin 2P + L 
2P 

where 

0 . „ h AA 

Pi tan p ± = 

Therefore, 


(A10) 


T(x,r ) = T n + ^ y~ P 

v * O I II.W 


sin p(l - x) 
h2 l=I sin 2p + 

2P 


r 

Jo 


G(A)e 


KPi f 

‘ h 2 


-(t-A) 


dA 


(x < h) 


(All) 


The present paper is concerned with a linear temperature-rise rate at the 
heated face and solutions can be obtained by substituting G(t) = Bt into 
equation (All) . However, by substituting the linear temperature-rise rate 
into equation (A8) and again finding an inverse transform by a formal 
application of the inversion integral, the following equation is obtained: 


T(x,t) 


= T 0 + B 


2h 2 K 


T + K — ^-(x - h) + — (x 2 - h 2 ) 
AA . 2 ' ' 


JPi 2 ^ 

P sin p(l - x) h 2 


3 


+ 


(A12) 


The series in this equation is rapidly convergent because of the appear- 
ance of in the denominator. 
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APPENDIX B 

ANALYTICAL APPROXIMATION WITH RADIATION AND CONIXJCTION 


A single procedure is presented in this appendix for analytically- 
approximating the temperature of the unheated face of a sandwich panel 
in which internal radiation is present. The following steps are involved: 

(a) Calculate the temperature distribution, considering conduction 

only. 

(b) Calculate the radiant heat transmission to the unheated face, 
assuming that the temperature distribution is that calculated in step (a). 
Convert this heat input into a temperature rise. 

(c) Correct the temperature of the unhe ated face for the radiant 
heat input. 

(d) Assume that the actual conduction to the unheated face is pro- 
portional to the difference between the temperature of the heated face 
and the temperature of the unheat'ed face as calculated in step ( c ) . 
Calculate the temperature rise of the unheated face caused by this heat 
input. 

(e) Add the temperature rises calculated In steps (b) and (d) to 
obtain the final approximation to the temperature of the unheated face 
of the panel. 

Corresponding to each of these steps are temperatures which will be 
designated as T a , Tb, T c , T^, and T g , respectively. The tempera- 
tures T a are those calculated by considering conduction only (eq. (All)). 
The radiant heat input to the unheated face ( step (b ) ) can be calculated 
by integrating equation (10), using the temperatures T a . This heat input 
can be converted to a temperature rise by dividing the integral of heat 
input by the heat capacity of the un heated face. This temperature rise 
is Tb . 

The temperature of the unheated face is now (step (c)): 

T c = T a + Tb 

From step (d), the heat conducted to the unheated face Is 

1 = j(Ti) a - (Tp) a " (®F)b] 
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and the temperature rise caused by conduction satisfies the relation 


q = cptp 


d (%)d 

dr 


or 


The temperature of the unheated face is the sum of the temperature rises 
due to conduction and radiation: 


T e = Tb + T d 

Figure 7 compares results calculated on the analog computer from equa- 
tion (lj) with results calculated hy the above procedure. The dashed 
line is the result obtained on the analog computer. The dash-dot line 
is (Tp) c , which has a correction for the effect of radiation but does 

not consider any interaction between radiation and conduction. The 
solid line is (Tp) e . The close agreement between the solid and dashed 
lines indicate that the approximate analysis is adequate up to the time 
at which the maximum temperature difference is reached. 
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TABLE I.- THE FIRST FIVE ROOTS OF THE EQUATION 0 ± tan 0 .J. 


h AA. 
tj> 


h AA 
tjr 

Pi 

P2 

^3 

P 4 

h 

O.Ol 

O.O998 

3.1448 

6.2848 

9.4258 

12.5672 

• 06 

.2425 

3.1606 

6.2927 

9 . 4311 

12.5711 

.1 

.3111 

3.1731 

6.2991 

9.4354 

12.5743 

.2 

.4328 

3.2039 

6.3148 

9.4459 

12.5823 

• 3 

.5218 

3.2341 

6.3305 

9.4565 

12.5902 

.4 

.5932 

3.2636 

6.3461 

9.4670 

12.5981 

.5 

.6533 

3.2923 

6.3616 

9.4775 

12.6060 

1.0 

.8603 

3.4256 

6.4373 

9.5293 

12.6453 


TABLE II.- OVERALL CONFIGURATION FACTORS 


Factor 

(a) 

1-0.8 

— — 1 2 

I - 1 - 6 

|- 2 .0 

ro 

H 

0.407 

0.527 

0.620 

0.676 

f 2,3 

.164 

.300 

.4oo 

.495 

f 2,4 

.092 

•133 

.145 

.125 

f 2,5 

.128 

.100 

.080 

.090 

F 3A 

.211 

.220 

.200 

.177 

F 5A 

.252 

.152 

.100 

.o4o 


Equivalent factors are: 


F 2 ,l = F 5,4 = F 4,5 
f 2,3 = F 3,2 = F 3,4 = Fj +,3 
f 2,4 = F 4,2 
f 2,5 = F 5,2 = F 4 ,l 
F 3 ,l = F 3>5 = F 5,3 


























STATION 2 
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o 


STATION 4 
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Time, sec 

Figure 2.- Typical calculated results, h = 0 .5 inch; tp = 0.03 inch; 
v = 15 Ib/ft3; B = 10° R per second. 
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clmum t emp erature difference with radiation and conduction 
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28 



Figure k.- Maximum dimensionless temperature difference as a function of 

h 


■F 
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Figure 4.- Continued. 
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Figure 5 .- Comparison of experimental and calculated results. 
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Figure 7- - Results of approximate analytical procedure. 
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